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Abstract 

We  describe  randomized  parallel  CREW  PRAM  algorithms  for  building  trapezoidal 
diagrams  of  line  segments  in  the  plane.  For  general  segments,  we  give  an  algorithm  re- 
quiring optimal  0(A  +  nlogn)  expected  work  and  optimal  O(logn)  time,  where  A  is  the 
number  of  intersecting  pairs  of  segments.  If  the  segments  form  a  simple  chain,  we  give 
an  algorithm  requiring  optimal  0(n)  expected  work  and  C)(logn  loglogn  log*  n)  expected 
time,  and  a  simpler  algorithm  requiring  0(n  log'  n)  expected  work.  The  serial  algorithm 
corresponding  to  the  latter  is  the  simplest  known  algorithm  requiring  0(nlog*  n)  expect- 
ed operations.  For  a  set  of  segments  forming  A'  chains,  we  give  an  algorithm  requiring 
0(A  +  nlog*  n  -I-  A'logn)  expected  work  and  0(  log  n  loglogn  log'  n)  expected  time.  The 
parallel  time  bounds  require  the  assumption  that  enough  processors  are  available,  with 
processor  allocations  every  logn  steps. 

1     Introduction 

1.1      Results  and  related  work 

We  give  several  algorithms  for  building  trapezoidaJ  diagrams  of  line  segments  in  the  plane. 
These  algorithms  reach  or  approach  optimality  with  respect  to  three  parameters:  the  number 
n  of  segments,  the  number  A  of  crossing  pairs  of  segments,  and  the  number  A'  of  chains.  Here 
a  chain  is  a  sequence  of  segments  such  that  a  segment  meets  its  successor  in  the  sequence  at 
an  endpoint,  and  no  endpoint  is  common  to  more  than  two  segments.  Specific  values  of  these 
parameters  correspond  to  problems  in  computational  geometry  that  have  been  extensively 
studied.  We  consider  these  cases,  as  well  as  the  general  case,  in  the  context  of  randomized 
parallel  algorithms  in  the  CREW  PRAM  model. 

Given  a  set  5  of  n  line  segments,  the  trapezoidal  diagram  T(S)  is  a  collection  of  simple 
regions  defined  as  follows:  for  each  endpoint  of  a  segment  in  5,  and  each  intersection  point 
of  segments  in  5,  extend  horizontaJ  segments  to  the  left  and  right  as  far  as  possible  without 
crossing  a  segment  of  5.  (Some  of  these  segments  will  be  rays,  in  fact.)  These  horizontal 
visibility  edges,  together  with  the  edges  of  S,  subdivide  the  plane  into  polygonal  cells  with 
four  or  fewer  sides,  that  in  general  are  trapezoids. 

Simple  polygons.  A  simple  polygon  is  bounded  by  a  set  of  noncrossing  edges  (line 
segments),  so  building  the  trapezoidal  diagram  for  a  simple  polygon  corresponds  to  the  special 
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case  where  ^4  =  0  and  A'  =  1.  Sometimes  a  triangulation  of  a  simple  polygon  is  desired, 
not  a  trapezoidal  diagram.  A  triangulation  of  a  simple  polygon  is  a  partition  of  its  interior 
into  triangles  whose  vertices  are  also  vertices  of  the  polygon.  Most  known  algorithms  for 
triangulation  begin  by  computing  a  trapezoidal  diagram  of  the  polygon's  boundary  segments. 
Given  the  trapezoidal  diagram,  a  triangulation  can  be  found  in  0(n)  operations,  both  serially 
[CI84,FM84]  and  in  O(logn)  time  in  parallel  [Goo89b]. 

Many  problems  involving  a  simple  polygon  can  be  solved  in  linear  time  if  a  triangulation  of 
the  polygon  is  known,  so  the  complexity  of  triangulation  is  a  particularly  appealing  question. 
Chazelle  recently  settled  this  question  by  finding  a  sophisticated  algorithm  requiring  Oin)  time 
in  the  worst  case  [Cha90].  An  earlier  Las  Vegas  algorithm  due  to  Clarkson,  Tarjan,  and  Van 
Wyk  requires  O(nlog*  n)'  expected  time  for  any  polygon  [CTVW89].  In  this  note  we  describe 
another  Las  Vegas  algorithm  requiring  expected  O(Tilog'n)  operations;  the  serial  version  is 
much  simpler  than  the  earlier  Las  Vegas  algorithm,  which  is  much  simpler  than  Chazelle's.  Sei- 
del  recently  discovered  a  similar  sequential  algorithm  [Sei].  A  parallel  version  of  our  algorithm 
requires  0(logn  loglognlog*  n)  expected  time,  with  0{n)  operations  if  Chazelle's  algorithm 
is  invoked  for  subproblems,  or  O(nlog'  n)  operations  if  not. 

Segments  with  many  chains.  In  the  case  where  A'  =  fi(n).  we  cannot  hope  to  do  better 
than  0(nlogn),  by  reduction  from  sorting.  If  we  assume  nondegenerate  input,  so  that  no  point 
meets  more  than  two  segments  except  at  endpoints,  the  parameter  A  is  also  the  output  size, 
so  we  can  do  no  better  than  Q(A  -\-  n\ogn)  operations.  This  has  been  achieved  in  the  serial 
case  both  deterministically  [CESS]  and  with  randomization  [Mul88,CSS9].  Recently  Goodrich 
[Goo89a]  has  found  a  parallel  algorithm  requiring  0(i4  +  n  log  n)  logn  operations  and  O(logn) 
time.  For  line  arrangements,  where  A  =  (j),  Hagerup  et  al.  [HJ\V90]  have  found  an  algorithm 
requiring  0(71^)  operations  and  O(logn)  expected  time.  If  two  sets  of  segments  B  and  C  are 
given,  with  the  segments  in  each  set  non-intersecting,  then  an  algorithm  of  Goodrich  et  al. 
[GSG89]  can  find  the  intersections  between  the  two  sets  with  0{A  +  n  log  n)  work  and  O(log  n) 
time,  where  n  is  the  total  number  of  segments  and  A  is  the  number  of  intersections.  Here  we 
use  these  results  to  show  that  OiA  +  n\ogn)  expected  operations  and  O(logn)  time  suffice  in 
general. 

The  general  case.  Recently  Chazelle  has  posed  the  problem  of  computing  the  trapezoidal 
diagram  of  a  non-simple  chain  of  segments,  that  is,  where  A'  =  1  but  A  is  not  necessarily  zero 
[Cha90].  We  generalize  this  question  to  allow  a  dependence  on  A'  >  1,  and  show  that  our 
techniques  permit  the  construction  of  T{S)  in  0(A  +  n  log*  n  -♦-  A'logn)  expected  work  and 
O(lognloglognlog*  n)  expected  time.  This  is  optimal  work  with  respect  to  A  and  A',  and 
within  log*  n  of  optimal  with  respect  to  n. 

1.2      Outline  of  the  paper 

The  remainder  of  the  introduction  has  a  subsection  devoted  to  technical  issues  regarding 
the  representation  of  the  output,  the  nondegeneracy  assumptions,  and  the  model  of  parallel 
computation  that  we  use.  There  follows  a  subsection  giving  the  probabilistic  lemmas  we  need 
in  the  analysis  of  our  randomized  divide-and-conquer  algorithms. 

Our  serial  algorithms  are  given  in  the  next  section;  the  parallel  algorithms  are  given  in  §3. 

'log*  n  is  the  least  t  such  that  log''*  n  <  1,  where  log*''  n  is  the  ith  iterate  of  the  logarithm:   log'"'  n  =  n, 
and  log'"'  n  =  loglog'""''  n  for  i  >  0. 


1.3      Some  technical  issues 

We  compute  representations  of  trapezoidal  diagrams  that  include  a  data  structure  for  navigat- 
ing between  adjacent  cells;  there  is  a  technical  issue  here  of  what  constitutes  adjacency.  We 
might  say  that  two  cells  are  adjacent  when  they  share  a  common  visibility  edge,  and  when  no 
input  segments  cross,  this  definition  is  adequate.  However,  when  input  segments  may  cross, 
we  must  also  be  able  to  navigate  between  any  two  cells  that  share  a  common  boundary.  We 
will  call  a  data  structure  recording  adjacency  under  the  former  definition  a  partial  adjacen- 
cy representation,  and  the  latter  a  complete  adjacency  representation.  A  complete  adjacency 
representation  might  include  for  each  cell  an  ordered  list  of  cells  incident  to  it. 

The  parameter  A  is  the  number  of  pairs  of  segments  that  cross:  that  is,  whose  nonempty 
intersection  is  not  an  endpoint  of  either  segment.  We  assume  in  general  that  no  point  meets 
three  segments  except  at  endpoints;  such  a  condition  is  easily  simulated  using  a  tie-breaking 
scheme. 

We  assume  each  chain  is  given  as  a  list  of  consecutive  line  segments,  in  clockwise  order, 
and  that  no  two  distinct  endpoints  have  the  same  y  coordinate;  the  latter  nondegeneracy 
assumption  is  also  easily  avoided  using  a  simple  tie-breaking  scheme.  Since  each  chciin  is 
readily  transferred  to  an  array  with  the  edges  in  chain  order,  we  assume  henceforth  that  the 
chains  are  stored  thus.  (In  the  parallel  case  this  transfer  is  done  by  a  list-ranking  algorithm 
[AM88];  it  requires  0(m)  operations  and  O(logm)  time  for  m  edges.) 

We  use  the  CREW  PRAM  model  for  parallel  computation,  with  random  integers  in  the 
range  1  to  n  available  at  unit  cost.  We  report  the  complexity  as  a  pair  {W,T),  where  W  is 
the  number  of  operations  (or  the  work)  done,  and  T  is  the  parallel  time  achievable  for  that 
number  of  operations.  This  contrasts  to  the  standard  form  of  results,  which  reports  the  pair 
{P,T),  where  P  is  the  number  of  processors  used  and  T  is  the  parallel  time.  Our  description 
eases  the  analysis  of  algorithms  where  the  number  of  processors  used  varies  over  the  course 
of  the  algorithm.  Of  course,  in  an  implementation,  the  number  of  processors  used  is  fixed; 
phases  that  require  more  processors  are  handled  by  a  round  robin  allocation  of  processors.  It 
has  to  be  verified  that  such  an  allocation  of  processors  is  feasible,  but  for  our  algorithms  this 
presents  no  difficulty. 

The  number  of  processors  (and  operations)  varies  for  two  reasons:  we  do  not  know  in 
advance  the  number  A  of  intersections  of  the  segments,  and  we  don't  know  the  results  of  some 
randomized  divide-and-conquer  steps.  Our  assumption  is  that  every  logn  steps  the  algorithm 
is  allowed  to  request  an  additional  allocation  of  processors,  which  will  then  be  provided  in 
part  or  all  by  the  system.  As  an  example  the  usefulness  of  this  assumption,  in  the  random- 
ized divide-and-conquer  steps,  we  obtain  a  collection  of  subproblems  with  total  expected  size 
within  appropriate  bounds;  we  know  the  subproblem  sizes,  and  allocate  processors,  before  the 
subproblems  are  solved. 

We  assume  further  that  the  processors  being  used  by  the  algorithm  are  always  numbered 
consecutively.  We  make  extensive  use  of  the  two  algorithms  mentioned  earlier,  due  to  Goodrich 
[Goo89a]  and  Goodrich,  Shauck  and  Guha  [GSG89].  While  neither  procedure  is  described  in 
exactly  this  format  it  is  not  difficult  to  rewrite  them  for  this  model.  Now  we  explain  the 
meaning  of  our  complexity  results.  The  parallel  time  assumes  that  processor  requests  are  met 
in  full;  the  algorithm  will  be  proportionately  slower  if  fewer  processors  are  provided;  regardless 
of  the  number  of  processors  provided  the  operation  counts  are  as  stated. 


1.4     Probabilistic  divide-and-conquer 

Our  aJgorithms  will  depend  on  the  following  properties  of  trapezoidal  diagrams  of  random 
subsets  of  S,  which  we  use  for  divide-and-conquer. 

Lemma  1  Given  a  set  S  of  n  line  segments  with  A  intersecting  pairs,  let  R  C  S  be  a  random 
subset  of  S  of  size  r,  with  all  subsets  of  size  r  equally  likely.  For  cell  T  E.  T{R),  let  uj  denote 
the  number  of  line  segments  of  S  that  meet  the  interior  of  T,  and  let  Nj  denote  the  number 
of  cells  adjacent  to  T  (sharing  a  boundary  point  or  edge).   Then: 

(i)  The  expected  number  of  intersections  of  R  is  no  more  than  Ar'^jn^. 

(ii)  The  expected  value  £'IIt6T(S)  ^tNt  =  0{n  -f  Arjn). 

(Hi)  The  expected  value  EYlT£T(S)'^T^ognT  =  0{(n  -\-  Ar/n)\og{n/r)). 

(iv)  With  probability  1  -  1/n,  max7-gx(S)  "T  =  0((logn)n/r). 

(v)  The  expected  value  EY,T£T[S)  "t  =  0{A  +  n'^/r). 

Proof  The  claim  (i)  is  Lemma  4.1  of  [CS89].  Claim  (ii)  can  be  proven  analogously  to 
Lemma  4.3  of  that  paper,  and  claim  (iv)  follows  from  Corollary  4.4  of  [CS89].  Claim  (iii)  is  a 
corollary  of  Theorem  3.6  of  that  paper,  with  c  =  2  and  VV{j)  -  y/j\ogj  in  the  notation  of  the 
theorem,  and  claim  (v)  also  follows  from  that  theorem  using  c  =  2  and  W{j)  =  _;'.    D 

2      Serial  Algorithms 

We  start  by  describing  an  algorithm  for  segments  that  form  a  simple  (non-self-intersecting) 
polygonal  chain,  then  the  case  of  a  single  polygonal  chain  that  may  have  self-intersections,  and 
finally  the  general  case  of  several  polygonal  chains. 

2.1      A  single  simple  chain 

The  algorithm  for  a  simple  polygonal  chain  is  a  modification  of  the  randomized  algorithm  of 
[CTVW89];  this  version  also  has  a  running  time  of  0{n  log"  n)  time;  its  novelty  is  that  it  does 
not  require  Jordan  sorting. 

Given  a  polygon  P  with  set  5  of  n  edges,  the  algorithm  builds  T{S)  in  log'  n  phases;  in 
phase  i,  the  diagram  T{S')  is  found,  where  5'  is  a  random  subset  of  S  of  size  r,  =  [n/[log*'^  n]\. 
In  fact  we  will  use 

5^  C  5^  C  •  •  •  C  S'*^*' "  =  5, 

with  5'  a  random  subset  of  5''*''  of  the  given  size;  imagine  the  segments  of  S  to  be  randomly 
permuted,  and  take  5'  as  the  first  r,  segments  in  the  permuted  order. 

Each  phase  will  require  an  expected  0(n)  operations,  as  discussed  below. 

The  diagram  T{S^)  can  be  computed  in  0{n)  time,  by  various  algorithms  that  require 
0(m  log m)  time  for  a  set  of  m  segments.  [CE88,Mul88,CS89]  This  completes  phase  1. 

Each  phase  i  >  1  has  three  parts: 

•  insert  the  segments  of  5  into  T{S'~^ );  that  is,  find  the  cells  that  meet  each  segment,  and 
so  the  set  of  segments  Sj  that  meet  each  cell  T  G  T(5'"')- 

•  for  each  cell  7  G  T(5'->),  compute  T  n  T{S'j),  where  5^  =  5t  n  S'. 


•  merge  appropriate  cells  of  the  diagrams  T  PI  T(Sj)  to  build  cells  of  T{S'). 

To  do  the  first  step  of  inserting  the  segments,  we  walk  along  the  segments  of  P,  and 
simultaneously  through  T{S'~^),  finding  for  each  segment  in  turn  the  cells  of  T{S^~^)  that 
it  meets.  By  the  nondegeneracy  assumption,  no  cell  shares  a  bounding  horizontal  edge  with 
more  than  four  others,  so  0(1)  time  is  required  for  each  cell  that  a  segment  meets. 

To  compute  the  diagrams  TnT(5j),  we  again  use  an  aJgorithm  requiring  0(m  log  m)  time 
for  a  set  of  m  segments. 

The  resulting  visibility  edges,  over  all  the  cells,  together  with  the  visibility  edges  for 
T{S'~^),  form  a  superset  of  such  edges  for  T(S').  Some  visibility  edges  of  T{S'~^)  are  not 
in  T{S'),  so  the  merging  step  of  phase  t  is  needed  to  obtain  the  cells  of  T(S').  Note  that  a 
group  of  cells  that  we  need  to  merge  are  in  some  collection  of  T  G  T(S'~^)  that  are  in  order 
along  a  segment  e  €  5',  and  we  know  this  order  because  it  was  obtained  during  the  walk  of 
the  polygon  and  T(S'~^). 

2.1.1      Analysis 

As  noted,  each  phase  of  the  algorithm  requires  0{n)  expected  time.  For  phase  i  >  1,  the 
expected  totaJ  time  of  the  insertion  step,  to  find  all  the  sets  5j,  is  proportional  to 

E      ^       \ST\  =  0(n) 

T6r(S'-') 

by  Lemma  l(ii).  Letting  nj  =  |5j|  (so  5'  plays  the  role  of  S  in  Lemma  1,  and  S'~^  the  role 
of  R),  the  expected  time  to  compute  T  n  T{S'j)  for  all  T  G  T{S'-^)  is 

E       Yl       "T  log  nr  =  0(15*1  log  ^1^) 

log^  '  n  log'  '  n 

=     0{n) 

by  Lemma  l(iii).  Finally,  the  merging  step  can  be  done  in  time  proportional  to  the  number 
of  trapezoids  in  all  the  T{Sj),  which  is  certainly  an  expected  0{n). 
We  conclude: 

Theorem  2  There  is  a  randomized  algorithm  that  triangulates  n-edge  simple  polygons  in  ex- 
pected 0{n  log'  n)  time. 

2.2      A  single  chain 

The  algorithm  requires  few  changes  when  the  chain  has  A  >  0  pairs  of  intersecting  segments, 
except  to  require  the  use  of  more  sophisticated  algorithms  to  compute  T{S^),  the  T{Sj  ), 
and  to  do  the  walks.  The  "base"  aJgorithm  for  trapezoidal  diagrams  must  take  0{k  +  Tn\ogm) 
expected  time  to  compute  the  trapezoidal  diagram  of  m  segments  with  k  intersecting  pairs. 

For  the  analysis,  it  suffices  to  bound  the  work  due  to  intersecting  pairs,  since  the  remaining 
work  is  0(n  log' n)  as  in  Theorem  2.  By  Lemma  l(i),  the  expected  number  of  intersections 
among  segments  of  5'  is  0{A/  log^  n),  which  bounds  the  expected  work  in  computing  T{S^). 


Since  the  work  in  traversing  a  cell  T  adjacent  to  TVj  other  cells  is  proportional  to  Nj,  the  work 
in  walking  the  chain  to  find  the  sets  5j  is  bounded,  by  Lemma  l(ii),  by  0{A/  \ogn)  (ignoring 
the  term  not  dependent  on  A).  In  general,  the  work  to  compute  the  intersection  points  of  S' 
is  expected  0(A/{\og^'^  ")^)>  which  is  dominated  by  the  time  to  traverse  T{S')  to  find  the  sets 
St,  which  is  0{A/ \og^'^  n).  The  expected  work  dependent  on  A  is  therefore 

l<.<log-nl°g"" 

which  is  0(A). 

2.3     Many  chains 

To  complete  this  section,  consider  the  general  case:  the  input  is  a  set  of  A'  chains  comprising  n 
segments,  having  A  intersecting  pairs.  We  modify  the  above  algorithm  only  slightly:  for  each 
chain,  choose  one  of  its  endpoints  to  be  the  leader  oi  the  chain.  The  algorithm  will  maintain 
the  cells  of  T(S')  contcuning  these  leaders,  and  use  planar  point  location  procedures  to  make 
this  operation  fast.  With  such  leaders  known,  the  walk  through  T{S^)  is  as  before,  starting 
with  the  leader  for  each  chain. 

Recall  that  a  planar  subdivision  defined  by  m  non-intersecting  straight  edges  can  be  pre- 
processed  in  O(mlogm)  time  so  that  the  region  containing  a  given  query  point  can  be  found  in 
O(logm)  time  (various  such  algorithms  are  discussed  in  [Ede87]).  We  apply  such  a  procedure 
to  T(S^),  and  find  the  cell  contciining  each  leader.  To  find  the  cell  of  T{S''^^ )  containing  each 
leader,  given  that  information  for  T(S'),  we  preprocess  each  TiSj"^)  for  point  location,  locate 
each  leader  within  the  appropriate  cell,  and  maintain  that  information  when  creating  T(5'"*'*). 

Theorem  3  Given  a  set  of  S  line  segments  forming  K  chains  and  with  A  intersecting  pairs 
of  segments,  the  trapezoidal  diagram  T{S)  can  be  found  in  0(A  +  nlog'  n  -f  A'logn)  expected 
time. 

Proof  Since  the  preprocessing  for  point  location  in  the  various  trapezoidal  diagrams  takes 
no  longer  than  the  construction  of  those  diagrams,  it  remains  only  to  find  the  time  needed  for 
the  queries. 

Locating  the  K  leaders  within  T{S^)  requires  C>( A'logn).  By  Lemma  l(iv),  for  i  >  0,  the 
maximum  size  of  a  set  5j"\  over  all  cells  T  6  T(S'),  is 

0(logn^^)  =  O(logn  log<'>  n/  log<'+')  n) 

with  probability  I  -  ^-  The  query  time  for  locating  leaders  within  T(S')  for  i  >  1  is  thus 
O(loglogn).  The  cost  of  all  queries  after  the  first  is  therefore  0(  A'loglogn  log*  n),  which  is 
dominated  by  the  cost  for  locating  within  T{S^).  If  some  trapezoid  T  G  T{S')  contains  more 
than  ©(log nlog*''  n/ log*'"*"'' n)  edges,  then  the  point  location  takes  at  most  O(logn)  time;  but 
this  happens  with  probability  at  most  ^,  and  so  makes  a  contribution  of  0(  A'logn  log' n/n) 
to  the  overall  expected  running  time,  a  negligible  term.    D 


3      Parallel  Algorithms 

We  describe  several  parallel  algorithms:  first,  one  suitable  for  A  =  n(n^),  then  for  K  =  i}{n) 
(or  at  least,  any  chain  connectivity  is  unused).  Following  this  is  an  algorithm  for  a  single 
simple  chain,  then  one  for  a  single  chain  with  intersections,  and  finally  the  general  case.  First 
we  discuss  the  use  of  random  subsets  for  divide-and-conquer,  in  the  parallel  setting. 

3.1  The  divide-and-conquer  scheme 

Each  parallel  algorithm  described  below  generally  has  one  or  more  phases  like  those  of  the 
serial  algorithm,  applying  randomization  for  divide-and-conquer: 

•  Compute  T{R)  for  random  R  C  S  of  size  r,  and  possibly  a  point  location  data  structure 
for  it; 

•  insert  the  segments  of  5  into  T(R);  that  is,  find  the  cells  that  meet  each  segment,  and 
so  the  set  of  segments  Sj  that  meet  each  cell  T  £  T(R)\ 

•  for  each  cell  T  £  T{R),  compute  T  fl  T{St)\  "  v  ■ 

•  merge  appropriate  cells  of  the  subdiagrams  T  DTiSj)  to  build  cells  of  T{S). 

The  first  step  generally  applies  some  previous,  possibly  suboptimal,  algorithm.  The  inser- 
tion step  is  more  problematic  for  parallel  algorithms  than  in  the  serial  case:  we  cannot  simply 
walk  along  chains.  In  fact,  single  segments  may  meet  ^(r)  trapezoids,  even  though  an  average 
segment  meets  no  more  than  6,  and  so  walking  through  T{R)  along  some  segments  is  too  slow. 
This  difficulty  is  the  main  problem  for  the  rest  of  this  paper. 

In  many  variations  of  the  insertion  step,  we  obtain,  in  effect,  a  collection  of  pairs,  comprising 
segments  of  S  and  cells  of  T{R)  that  intersect.  We  do  not  have,  for  each  cell  T  E  T(R),  a.  list 
St  of  segments  of  5  that  meet  it.  To  create  such  lists,  the  pairs  are  integer  sorted,  using  their 
trapezoid  as  the  key.  Such  sorting  can  be  done  with  the  algorithm  of  Rajasekaran  and  Reif 
[RR88]  in  O(logn)  time  and  work  linear  in  the  number  of  sorted  pairs. 

As  in  the  sequential  aJgorithms,  the  merge  step  builds  a  trapezoid  T  of  the  finaJ  diagram 
from  several  trapezoids  in  the  subdiagrams.  This  is  done  as  follows.  Implicitly,  for  each  edge 
of  S,  the  trapezoidal  diagrams  (for  each  cell  T  6  T{R))  define  a  linked  list  of  its  intersections 
in  order.  By  applying  a  list  ranking  algorithm  to  this  list  [AM88,CV86]  we  eliminate  the 
intersections  with  visibility  edges  of  T(R).  This  leaves  the  intersections  among  the  edges 
of  5  together  with  the  visibility  edges  of  T{S);  the  computation  of  T{S)  is  completed  by 
linking  adjacent  edges.  The  complete  adjacency  representation,  which  by  Lemma  l(ii)  is 
of  size  0{n  -\-  A),  is  obtained  by  a  further  application  of  list  ranking  to  the  edges  of  S  so 
as  to  determine  for  each  trapezoid  the  number  of  its  neighbors;  a  proportionate  number  of 
processors  is  then  allocated  to  each  trapezoid;  these  processors  are  used  to  create  the  adjacency 
lists,  each  processor  being  responsible  for  copying  the  name  of  one  trapezoid  to  the  adjacency 
list.  Overall,  this  takes  0{n  -\-  A)  operations  and  O(logn)  time. 

3.2  General  line  segments 

This  subsection  gives  an  algorithm  supporting  the  following  theorem. 


Theorem  4  Given  a  set  S  of  n  line  segments  in  the  plane  with  A  intersecting  pairs,  the 
trapezoidal  diagram  T(S)  can  be  found  in  0{A  +  nlogn)  expected  operations  and  in  O(logn) 
worst-case  time  on  a  CREW  PRAM. 

The  algorithm  uses  the  divide-and-conquer  approach  described  above:  first,  find  an  esti- 
mate A  o{  A.  Then  take  a  random  subset  Rc  S  of  size  r  =  n^/{A  +  nlogn).  Compute  T{R) 
and  process  it  for  point  location  using  previous  algorithms  [Goo89a,ACG89];  then  for  each 
cell  T  €  T(R),  find  the  sets  St  of  segments  that  meet  it,  and  compute  visibility  information 
among  the  segments  of  Sj  using  an  adaptation  of  the  algorithm  of  Hagerup  et  al.  [HJVV90]. 
Finally,  merge  the  resulting  cells  to  form  T{S). 

The  idea  is  to  balance  between  these  known  algorithms,  building  the  relatively  small  T{R) 
to  give  subproblems  with  many  intersecting  pjiirs  relative  to  the  number  of  segments,  so  that 
an  algorithm  like  that  of  [HJW90]  is  economical. 

Here  is  the  algorithm  in  more  detail:  the  estimate  A  is  found  by  taking  n  random  pairs 
of  segments  and  counting  the  number  of  intersecting  pairs  among  them.  Each  pair  (a,,  6,)  is 
chosen  from  among  the  (j)  pairs  of  segments  of  5,  with  all  pairs  equally  likely.  If  this  count 
is  C,  take  A  =  (C  +  l)(2)/n.  (Other  methods  for  estimating  A  might  well  be  acceptable,  but 
this  is  easy  to  analyze.)  In  particular,  we  note  that 

Lemma  5    The  expected  value  EA  =  A  +  (j)/",  and  E[l/ A]  <  l/A. 

Proof.  The  quantity  C  is  a  binomial  random  variable  with  n  trials  and  success  rate  per 
trial  p=  /l/(5).  We  have 

and  since  (p(n  +  l)/(j  +  1)  =  Q+J),  we  have 

which  is  no  more  than  1/pn,  or  {'^)/An,  and  so  ^^[l/^]  <  l/A.    D 

Insertion.  The  insertion  step  of  finding  the  sets  Sj  for  cells  T  £  T{R)  is  done  as  follows: 
process  T{R)  for  planar  point  location  [ACG89],  locate  the  endpoints  oi  S  \  R  within  T{R), 
and  then  in  parallel,  walk  along  each  segment  through  the  diagram,  finding  cells  that  it  meets. 
This  is  done  until  all  but  n/logn  segments  are  completely  traversed:  do  the  walk  log  n  steps, 
check  the  number  of  segments  not  yet  completed,  walk  along  logn  steps,  check  again,  and  so 
on,  until  the  condition  holds. 

There  remain  at  most  n/logn  segments  not  completely  inserted,  that  is,  not  completely 
traversed.  To  insert  these,  we  split  them,  at  random,  into  logn  groups  each  of  size  at  most 
n/  log^  n.  The  algorithm  of  [Goo89a]  is  applied  to  each  group  to  find  all  its  intersection  points, 
and  to  split  the  segments  of  the  group  up  into  segments  that  meet  only  at  endpoints.  For  each 
group,  we  find  the  intersections  of  the  resulting  segments  with  the  visibility  edges  of  T{R) 
using  the  algorithm  of  [GSG89]. 

We  now  have  the  pairs  of  segments  of  5  and  trapezoids  of  T{R)  that  meet;  as  discussed 
in  §3.1,  we  now  use  integer  sorting  to  obtain  the  list  St  of  segments  that  meet  a  trapezoid  T, 
for  every  T  £T{R). 
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Subdiagrams.  Next,  for  each  cell  T  G  T{R),  we  compute  the  subdiagrams  T  n  T{St) 
using  an  adaption  of  the  algorithm  of  [HJW90].  Our  adaption  takes  0(ti'^)  expected  opera- 
tions to  compute  the  trapezoidal  diagram  of  a  set  of  n  line  segments,  and  requires  O(logn) 
time.  While  their  algorithm  finds  triangulations  of  arrangements  of  lines,  we  are  interested 
in  trapezoid^d  diagrams  of  line  segments.  Nonetheless,  their  algorithm  can  be  adapted  to  our 
purpose;  Appendix  H  gives  the  deteiils. 

The  merge  step  proceeds  as  discussed  in  §3.1. 

Proof  of  Theorem  4.  The  proof  of  Theorem  4  is  completed  by  analyzing  the  above  al- 
gorithm. By  Lemma  l(ii),  the  expected  work  in  walking  along  the  segments  through  T{R) 
is  0{n  +  An/{A  +  nlogn)),  excluding  busy  waiting  by  processors  handling  segments  whose 
traversal  is  complete.  But  E[A/A]  <  1,  so  the  expected  work  is  0{n).  This  implies  that  after 
some  clogn  steps,  for  a  large  enough  constant  c,  at  most  n/  logn  segments  are  not  completely 
traversed. 

The  cost  of  inserting  these  remaining  segments  is  bounded  as  follows.  The  algorithm  of 
[Goo89b]  requires  an  expected  0{n\ogn  +  A)  operations  and  O(logn)  time  (for  by  Lemma  l(i), 
there  are  an  expected  0{A/  logn)  intersections  overall  present  in  these  randomly  selected 
subsets),  while  the  algorithm  of  [GSG89]  requires  an  expected  0( nlogn  +  -4)  operations  and 
O(logn)  time  (the  expectation  arises  because  there  are  an  expected  0{n  +  A/ logn)  edges 
present  over  all  the  subproblems). 

The  expected  number  of  operations  for  all  executions  of  the  algorithm  of  Hagerup  et  al. 
is  0(IIr€T(fi) '^r)!  ^y  Lemma  l(v),  this  has  expected  value  0{A  +  n^/r)  with  respect  to  the 
choice  of  R,  and  the  expected  value  of  n^ /r  with  respect  to  the  estimation  of  A  is  0(A  +  n  log  n), 
using  Lemma  5.  The  running  time  of  this  step  is  O(logn). 

The  complexity  of  the  final  step  (merging  trapezoids)  is  bounded  as  follows.  By  Lem- 
ma l(ii),  there  are  an  expected  0{n  +  Ar/n)  =  0(n)  intersection  points  along  the  edges  of 
S  (including  intersections  with  the  visibility  edges  of  T{R)).  So  the  complexity  of  the  prefix 
lists  algorithm  is  an  expected  0(n)  operations  and  O(logn)  time. 

On  summing  the  complexities  of  each  step,  the  theorem  follows  readily.    D 

3.3     A  single  simple  chain 

The  algorithm  proceeds  in  the  same  log*  n  phases  as  the  sequential  algorithm;  each  phase  uses 
the  randomized  divide-and-conquer  scheme,  with  T{S'~^)  used  to  build  T{S')  in  phase  i  >  1. 

The  diagram  T{S^)  is  built  in  phase  1  using  the  algorithm  of  Atallah,  Cole  and  Goodrich 
[ACG89];  it  does  0{n)  operations  in  time  O(logn).  Also,  we  create  a  planar  point  location 
structure  for  T{S^),  again  using  an  algorithm  of  Atallah  et  al. 

First  we  consider  phase  2,  where  we  build  T{S^)  in  parallel  using  an  expected  0(n)  oper- 
ations and  O(logn  log  logn)  time. 

Insertion.  The  insertion  of  edges  for  phase  2  is  done  by  traversal  of  T(5^)  and  the  edges  of 
the  chain,  as  in  the  sequential  algorithm.  However,  we  do  many  such  traversals  independently 
from  several  places. 

Initially,  the  chain  is  partitioned  into  subchains  of  length  logn.  This  takes  0(n)  operations 
and  O(logn)  time.  There  are  then  at  most  2n/logn  subchains.  We  use  n/(logn  log  logn) 
processors.  We  proceed  in  stages,  continuing  until  at  most  n/logn  edges  have  yet  to  be 
completely  inserted.  A  stage  lasts  logn  steps.  Each  step  comprises  the  traversal  of  an  edge  to 
its  next  intersection  point  or  to  its  far  endpoint,  whichever  is  nearer.  At  the  end  of  a  stage, 
if  fewer  than  ^n/(logn  log  logn)  chains  remain,  the  chains  are  evenly  partitioned  to  obtain 


between  ^n/(lognloglogn)  and  n/(lognloglogn)  chains  anew  (since  the  chains  are  contiguous 
sequences  of  edges  of  the  polygon,  storing  these  edges  in  clockwise  order  in  an  array  makes 
the  partitioning  straightforward).  The  new  endpoints  are  located  using  the  point  location 
data  structure,  at  a  cost  of  O(logn)  operations  per  endpoint.  So  a  stage  does  O(n/loglogn) 
operations.  It  either  detects  at  least  \nl  loglog  n  intersection  points  and  endpoints,  or  it  halves 
the  size  of  each  remaining  chain  (or  possibly  both).  By  Lemma  1,  there  are  an  expected  0{n) 
intersection  points,  so  after  an  expected  O(loglogn)  stages,  each  chain  comprises  at  most  one 
edge,  so  there  are  at  most  n/  logn  edges  which  have  not  been  completely  inserted;  that  is,  the 
procedure  terminates  within  an  expected  O(loglogn)  stages. 

To  insert  these  remaining  bad  edges,  we  find  the  intersections  between  them  and  the 
2n/logn  visibility  edges.  The  intersections  among  these  segments  are  computed  using  the 
algorithm  of  Goodrich  et  al.  [GSG89];  for  our  application,  it  requires  an  expected  0(n)  op- 
erations and  O(logn)  time  (the  expectation  arises  because  the  number  of  intersections  is  an 
expected  0{n)). 

As  discussed  in  §3.1,  we  now  use  integer  sorting  to  collect  the  sets  St  for  each  T  G  T{S^), 
using  the  algorithm  of  [RR88]. 

Subdiagrams.  Next,  for  each  T  G  T{S^),  we  find  T  n  T{Sl),  and  also  build  a  point 
location  data  structure  for  it,  with  the  procedures  of  Atallah  et  al.  [ACG89].  (The  point 
location  structure  will  be  used  for  later  processing.)  We  use  |5j-|  processors  for  this  task. 

Merging.  Finally,  we  need  to  merge  cells  from  among  those  of  the  subdiagrams  to  find 
those  of  T(S^).  To  carry  out  this  process,  for  each  edge  in  5^,  it  is  useful  to  have  a  linked  list 
of  its  intersection  points  with  the  visibility  edges  of  T(5' ).  Additionally,  in  phase  3,  it  will  be 
useful,  for  each  edge  in  5,  to  have  in  sorted  order  its  intersections  with  the  visibility  edges  of 
T{S^).  So  next  we  compute  for  each  edge  in  5  a  list  of  its  intersections  with  T(5'),  in  sorted 
order  along  the  edge. 

It  suffices  to  sort  the  intersection  points  for  the  0{nl  logn)  bad  edges,  since  the  other  inter- 
section points  are  already  in  order  aJong  the  edges  on  which  they  lie.  First,  we  randomly  select 
n/logn  of  these  intersection  points.  For  each  bad  edge,  we  sort  the  selected  points  along  the 
edge  and  split  the  edge  into  edge  portions  accordingly  (this  requires  two  sorts);  sorting  takes 
0{n)  operations  and  O(logn)  time  [Col88].  Next,  we  traverse  the  edge  portions,  seeking  inter- 
section points;  each  edge  portion  is  traversed  until  either  -logn  log  logn  intersection  points 
are  found,  or  the  portion  is  completely  traversed,  whichever  occurs  sooner;  a  is  a  constant  to  be 
specified.  The  edge  portions  are  evenly  redistributed  among  the  n/(logn loglog n)  processors 
every  logn  traversal  steps,  called  a  stage.  For  any  stage,  either  at  least  ^n/ log  logn  intersec- 
tion points  are  traversed,  or  the  processing  of  at  least  one  half  of  the  remaining  edge  portions 
is  completed.  Now  we  analyze  this  traversal:  the  expected  number  of  edge  portions  containing 
more  than  clognloglogn  intersection  points  is  (9(n/(logn)''''®*'^')  =  O(n/(logn)'"'''*'^),  say. 
It  follows  that  the  edge  portions  containing  more  than  -^  logn  log  logn  intersection  points  be- 
tween them  contain  an  expected  O(n/logn)  intersection  points.  So  the  procedure  lasts  an 
expected  C>(loglogn)  stages,  which  is  O(logn  log  logn)  time;  an  expected  0{n)  operations  are 
done.  In  addition,  there  now  remain  an  expected  0( n/logn)  untraversed  intersection  points. 
These  points  have  already  been  determined  in  the  previous  paragraph;  they  can  now  be  sort- 
ed with  respect  to  the  bad  edges  on  which  they  lie  and  then  merged  with  the  intersection 
points  Jilready  on  these  edges.  Overall,  sorting  the  intersection  points  takes  an  expected  0{n) 
operations  and  expected  O(logn  log  logn)  time. 

With  this  information,  the  merging  step  can  be  done  as  discussed  in  §3.1. 

The  algorithm  for  finding  the  diagrams  in  phases  t  >  2  is  analogous.   Before  we  describe 
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this  step,  here  is  analysis  of  resource  bounds  for  building  T{S^):  as  noted,  the  insertion  step 
requires  O(lognloglogn)  expected  time  and  0(n)  expected  operations,  and  building  the  lists 
57-  requires  no  more  resources.  The  construction  of  the  diagrams  TnT{Sj)  for  all  T  G  T{S^) 
requires  expected 

0{     Yl     \Sh^og\Sl\)  =  0{\S'\\og^-^)  =  0{n) 
rer(si)  '     ' 

operations,  and  O(logn)  time.  Sorting  the  intersection  points,  as  already  observed,  requires 
an  expected  0(n)  operations  and  O(lognloglogn)  time.  The  list  ranking  algorithms  require 
0(m)  operations  and  O(logm)  time  for  a  list  of  length  m,  so  the  expected  number  of  operations 
to  merge  cells  is  only 

0{     Yi      \S'T\)  =  0{\S'\)  =  0{n/\og\ogn). 

We  turn  to  the  procedure  for  building  T[S^)  given  T{S^~^),  for  i  >  2.  The  first  step  is  to 
insert  the  edges  of  the  polygon  into  T(S'~^).  The  edges  are  partitioned  at  the  points  at  which 
they  cross  visibility  edges  of  T(5'~^);  this  increases  the  number  of  edges  to  an  expected  0{n). 
The  same  edge  insertion  method  is  used  as  for  the  case  i  =  2.  Following  the  partitioning  of 
chains  into  subchains  of  length  logn  there  are  an  expected  0(n/  log''~^'  n)  chains  at  hand.  As 
before,  we  proceed  in  stages  until  at  most  n/logra  edges  have  not  been  completely  inserted. 
The  only  detail  of  the  insertion  that  is  novel  is  that  for  each  endpoint  of  each  edge,  the 
cell  T  6  T{S'~^)  containing  the  endpoint  is  known;  thus  a  point  location  on  an  endpoint 
requires  a  search  in  the  point  location  data  structure  for  T.  Note  that  only  0(n/logn)  point 
locations  are  done,  as  before.  The  insertion  of  the  n/  log  n  bad  edges  proceeds  essentially  as 
before,  except  that  a  separate  subproblem  is  created  for  each  cell  of  T{S'~^).  However,  the 
processor  allocation  and  reallocation  is  done  globally  across  all  these  subproblems;  this  allows 
for  variations  in  the  number  of  intersections  found  in  the  various  subproblems  and  allows 
expectations  to  be  with  respect  to  the  sum  of  the  sizes  of  the  subproblems.  Then,  for  each  cell 
T  G  T{S'~^),  we  obtain  a  list  of  the  segments  S  that  intersect  it,  as  before. 

Next,  for  each  cell  T  G  T{S'-^),  we  find  T  f)  T(S'j)  and  also  build  a  point  location  data 
structure  for  it,  using  the  procedure  of  Atallah  et  al.  [ACG89].  We  use  \Sj\  processors  for  this 
task. 

Finally,  the  cells  are  merged  and  for  each  edge  the  order  of  its  intersection  points  with 
T{S'~^)  is  computed  as  before. 

The  analysis  is  similar  to  that  for  the  case  i  =  2.  Again,  the  insertion  step  requires 
O(logTiloglogn)  expected  time  and  Oin)  expected  operations.  The  only  difficult  part  is  to 
analyze  the  cost  of  inserting  the  bad  edges.  Let  n'j  denote  the  number  of  bad  edges  intersecting 
cell  T  G  T{S'~^).  Then  the  operation  count  for  inserting  the  bad  edges,  using  the  algorithm 
of  Goodrich  et  al.  [GSG89],  is 

0(    E    {\S'T-'\  +  n'T)\og{\S'f'\  +  n'j)) 
=    0(|5-'|log||^+     Yl    "rlognV) 

This  is  an  expected  0(n),  since  '^t€T{S'-^)'"''t  -  0{n/\ogn).  Building  the  lists  5j  requires 
an  expected  0{n)  operations  and  O(logn)  time.    The  construction  of  the  diagrams  T(Sj-) 
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requires  expected 

0{      ^       |5^|log|5i^|)  =  0(|5'|)Iog  J^  =  0(n) 

operations  and  O(logn)  time.  Sorting  the  intersection  points,  as  before,  requires  an  expected 
0{n)  operations  and  O(lognloglogn)  time.  The  expected  number  of  operations  for  merging 
cells  is  only 

0{      Y.       |5^l)  =  0(l'?'l)  =  0(n/log<-)n). 

Summing  over  all  phases,  we  obtain  that  the  overall  expected  number  of  operations  done  is 
0(n  log*  n)  and  the  expected  time  required  is  O(lognloglog  n  log*  n).  This  constructs  the  hor- 
izontal visibility  structure.  It  remains  to  obtain  the  triangulation  from  the  visibility  structure. 
But  this  can  be  done  in  O(logn)  time  and  0{n)  operations  ([Goo89b]).  We  have  shown 

Theorem  6  There  is  a  parallel  randomized  CREW  PRAM  algorithm  for  computing  a  trian- 
gulation of  a  simple  polygon  which  runs  in  expected  0( log  n  log  log  n  log*  n)  time  and  does  an 
expected  0(n  log*  n)  operations. 

A  more  complex  algorithm  which  does  an  expected  0(n)  operations  in  expected  O(lognloglogn) 
time  can  be  obtained;  it  uses  Chazelle's  linear-time  triangulation  algorithm  as  a  subroutine. 

Here,  T{S^)  is  constructed  as  before  and  the  edges  of  the  polygon  are  inserted  into  T(S'^). 
By  Lemma  l(iv),  with  probability  at  least  1  -  ^,  each  cell  T  G  T(5^)  contains  at  most 
log  n  log  log  n  edges.  Using  the  method  of  Clarkson  et  al.  [CTV\V89]  the  edges  crossing  each 
visibility  segment  of  T(5'^)  are  Jordan  sorted  by  a  sequential  algorithm  (this  can  be  done  using 
Chazelle's  observation  that  Jordan  sorting  can  be  reduced  in  linear  time  to  triangulation,  and 
applying  Chazelle's  linear-time  triangulation  algorithm);  this  then  creates  a  set  of  visibility 
subproblems  which  are  solved  sequentially  by  applying  Chazelle's  linear-time  algorithm.  S- 
ince  each  subproblem  is  of  size  O(logTiloglogn)  with  overall  probability  1  -  i,  this  takes  an 
expected  0{n)  operations  and  expected  O(lognloglogn)  time.  This  then  yields 

Theorem  7  There  is  a  parallel  randomized  CREW  PRAM  algorithm  for  computing  a  trian- 
gulation of  a  simple  polygon  which  runs  in  eipec<edC)(lognloglogn)  time  and  does  an  expected 
0{n)  operations. 

3.4     A  single  chain  with  intersections 

The  algorithm  for  a  single  chain  with  A  >  0  intersections  is  similar  to  that  for  a  chain  with 
no  intersections. 

We  begin  by  computing  an  estimate  A  for  A  as  in  the  algorithm  of  §3.2.  If  A  >  nlogn, 
then  we  use  the  algorithm  given  in  Theorem  4  to  compute  the  trapezoidal  diagram  induced 
by  the  chain. 

Otherwise,  we  define  sets  5',5^,  •  •  • ,  5'°*'"  as  in  the  algorithm  for  a  simple  chain.  We 
begin  by  determining  T{S^),  but  now  we  use  the  randomized  algorithm  given  in  §3.2.  By 
Lemma  l(i),  the  expected  number  of  intersections  among  segments  of  5'  is  0{A/\og^  n), 
so  the  randomized  algorithm  of  §3.2  requires  an  expected  OiA/log^n  +  n)  operations  and 
O(logn)  time.  Then  we  insert  the  edges  of  the  chain  in  T{S^)  as  in  §3.3.  The  one  change  to 
the  traversal,  as  in  the  sequential  case,  is  that  when  and  if  the  traversal  crosses  a  segment  of 
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S^  the  work  required  is  proportional  to  the  number  of  cells  adjacent  to  this  edge  of  the  cell. 
So  to  analyze  the  traversal,  we  count  both  the  number  of  intersections  found  and  the  number 
of  cells  considered  in  crossing  segments  of  5';  by  Lemma  l(ii),  this  totals  an  expected  0(n), 
and  so  after  an  expected  O(loglogn)  iterations,  at  most  n/logn  edges  are  not  fully  inserted. 
A  new  method  is  needed  to  find  the  intersection  points  involving  these  bad  edges,  for  it  is  not 
guaranteed  that  the  bad  edges  do  not  intersect.  Instead  of  using  the  algorithm  of  Goodrich  et 
al.  [GSG89],  the  randomized  parallel  algorithm  given  in  §3.2  is  used.  To  apply  it,  as  in  §3.2, 
the  n/  log  n  bad  segments  are  randomly  partitioned  into  log  n  groups  of  n/  log  n  segments.  By 
Lemma  l(i),  each  group  of  segments  has  an  expected  0(  A/ log^  n)  intersections.  The  algorithm 
of  §3.2  is  applied  separately  to  each  group  of  segments  with  the  edges  of  T{S^).  Over  all  logn 
groups  the  algorithm  requires  an  expected  0(n  +  A/ logn)  operations  and  O(logn)  time  (for 
there  are  an  expected  0{n  +  /I/ logn)  intersection  points  over  all  logn  subproblems). 

The  sorting  of  intersection  points  is  done  as  in  §3.3. 

For  each  cell  T  £  T{S^)  we  obtain  a  list  of  the  segments  intersecting  T  as  in  §3.1.  Overall 
this  first  phase  requires  an  expected  0{n  +  A/logn)  operations  and  expected  O(lognloglogn) 
time. 

The  algorithm  proceeds  as  in  the  simple  chain  case,  as  modified  four  paragraphs  above, 
until  n/log'-'+^'n  =  |5-'+'|  >  nV(n  +  M  (i.e.  nlog'-''n  >  n  +  A  >  nlog<-'+'' n).  So  for  each 
cell  T  e  T{S^)  the  set  St  of  segments  intersecting  T  has  been  computed.  For  if  we  tried 
to  proceed  further,  the  insertion  of  the  chain  in  T(5-''*"')  would  be  too  expensive.  Then  the 
randomized  algorithm  of  §3.2  is  applied  to  each  set  St-  The  expected  work  for  this  final 
phase  is  0{YiTeT(SJ)  njlognj  +  A)  which  by  Lemma  l(iii)  is  0{n  +  |5-'|/l/n)log(n/|5-'|)  = 
0(n  log'-'"''''  n)  =  0{n  +  A);  the  running  time  is  O(logn). 

Phases  1  through  j  each  have  an  expected  operation  count  of  0{n  +  A/\ogn)  and  an 
expected  running  time  of  O(lognloglogn). 

Overall,  we  obtain 

Theorem  8  There  is  a  parallel  randomized  CREW  PRAM  algorithm  for  computing  the  trape- 
zoidal diagram  of  a  chain  with  A  intersections  which  runs  in  expected  0(\ogn\og\ognl3{  A,  n)) 
time  and  does  an  expected  0{A  +  n/3(/l,n))  operations,  where  l3(A,n)  =  minj  such  that 
A  >  n  log'-''  n,  for  A  >  n,  and  l3(A,  n)  —  log*  n  for  A  <  n. 

3.5      Many  chains 

As  in  the  serial  case,  consider  a  set  of  n  segments  comprising  K  chains,  having  A  intersecting 
pairs.  For  each  chain,  choose  one  of  its  endpoints  to  be  the  leader  of  the  chain.  The  algorithm 
maintains  the  cells  of  T{S')  containing  these  leaders,  using  planar  point  location  procedures 
to  make  this  operation  fast.  The  walk  through  T{S')  proceeds  as  before,  starting  at  the  leader 
of  each  chain.  Since  a  planar  point  location  data  structure  has  already  been  computed  for 
T{S^),  and  for  T(5'"'"')  with  respect  to  each  trapezoid  of  T(5'),  the  additional  expected  work 
for  doing  the  point  locations  on  the  leaders  is  just  C>(A'logn),  as  in  §2.3,  and  takes  expected 
time  O(logn).  (As  in  §2.3,  with  probability  at  most  ^,  a  point  location  takes  time  O(logn), 
and  otherwise  takes  time  O(loglogn)  for  :'  >  1;  overall,  this  is  an  expected  O(Iogn).)  So  we 
obtain: 

Theorem  9  There  is  a  parallel  randomized  CREW  PRAM  algorithm  for  computing  the  trape- 
zoidal diagram  of  K  chains  with  A  intersections  which  runs  in  expected 0{\ogn  \og\ogn/3{A,  n)) 
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time  and  does  an  expected  0{A  +  n0(A,n)  +  A'logn)  operations,  where  ^{A,n)  =  min;  such 
that  A  >  n  log^-'^  n,  for  A  >  n,  and  f3{A,  n)  =  log'  n  for  A  <  n. 
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Appendix  H 

We  have  a  set  5  of  n  line  segments,  and  we  want  to  find  the  trapezoidal  diagram  of  S  in  0{n^) 
expected  time.  The  algorithm  uses  randomized  divide-and-conquer:  take  random  R  C  S  oi 
size  r  =  n/logn,  and  compute  T(R)  using  Goodrich's  algorithm  [Goo89a].  (Alternatively, 
use  an  algorithm  anaJogous  to  the  suboptimaJ  one  given  in  [HJW90].)  Insert  the  remaining 
segments  of  S  into  T(R),  and  use  subsampling  of  the  segments  meeting  each  T  £  T{R)  to  yield 
subproblems  all  of  size  0(\/log  n)  with  high  probability.  Now  use  an  optimal  serial  aJgorithm 
for  such  subproblems,  and  Goodrich's  algorithm  for  the  remainder. 

The  insertion  step,  in  more  detail:  note  that  the  complete  adjacency  representation  of 
T{R)  divides  the  segments  of  R  into  0{r'^)  total  pieces,  with  O(r^)  on  any  segment  of  R 
(actually  0(ra{r))  trapezoids  meet  any  given  segment,  but  the  larger  bound  will  suffice).  We 
use  list  ranking  [AM88]  to  obtain,  for  each  segment  in  R,  a  sorted  array  of  the  pieces  induced 
by  T{R),  with  pointers  to  the  trapezoids  meeting  each  piece.  Now  for  each  segment  a  £  S, 
use  binary  search  on  the  array  for  every  segment  b  E  R,  to  find  the  trapezoids  containing  the 
intersection  point  of  a  and  b. 

This  does  not  give  all  the  trapezoids  meeting  a;  the  remaining  trapezoids  are  obtained 
as  follows:  consider  the  convex  diagram  of  R,  the  subdivision  induced  using  visibility  edges 
from  segment  endpoints  only,  not  intersection  points.  The  trapezoids  within  each  region  of  the 
convex  diagram  can  be  ordered  top  to  bottom;  such  orderings  can  be  obtziined  by  list  ranking, 
applied  to  adjacencies  between  trapezoids  with  common  edges  that  are  visibility  segments 
from  intersection  points.  The  trapezoids  in  a  convex  cell  that  meet  a  segment  a  G  5  are  an 
interval  in  that  list,  and  the  highest  and  lowest  trapezoids  contain  either  intersection  points 
of  a  with  R,  or  endpoints  of  a.  Binary  searches  suffice  to  locate  the  endpoints  of  a  in  the  list, 
if  either  are  contained  in  trapezoids  in  the  cell.  Finally,  the  trapezoids  meeting  a,  other  than 
the  highest  and  lowest,  can  be  obtained  after  sufficient  processors  are  allocated  to  do  this. 
(Note  that  the  number  of  trapezoids  met  is  available.)  We  have,  for  every  a  E  S,  the  set  of 
trapezoids  of  T(R)  that  it  meets.  By  integer  sorting  we  obtain  a  collection  of  subproblems  as 
discussed  in  §3.1.  This  completes  the  first  phase  of  processing. 

Now  for  each  T  €  T(R),  take  a  random  subset  of  the  segments  that  meet  it,  of  size 
A'|n7-|  log  |n7-|/v^ogn,  where  nj  is  the  number  of  segments  meeting  T,  and  K  is  an  appropriate 
constant.  Agjiin  apply  Goodrich's  algorithm  to  each  such  subset,  and  use  the  same  insertion 
technique,  to  obtain  a  collection  of  subproblems.  Now  after  appropriate  processor  allocations, 
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apply  an  optimal  serial  algorithm  to  those  subproblems  with  no  more  than  y/\ogn  segments, 
and  apply  Goodrich's  algorithm  to  the  remainder. 

Analysis.  It  is  easy  to  verify  that  the  algorithm  consists  of  a  constant  number  of  stages 
each  requiring  O(logn)  time  in  the  worst  ca^e. 

It  is  also  easy  to  verify  that  0{n^)  expected  work  is  required,  using  Lemma  1.  For  example, 
computing  the  trapezoidal  diagrams  of  the  random  subsets  in  the  second  phase  of  processing 
requires  expected  work  proportional  to 

^    [0(|nr|)log|nr|/yiogn]Mog|nr|; 

TeT(R) 
since  |n7-|  =  O(log^  n)  for  all  T  with  probability  1  -  1/n,  the  above  is  no  more  than 


J2     0(|nT|')(loglogn)Vlogn, 
TgT{R) 


which  is  0(n^(loglogn)^/log7i). 
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